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Abstract 

We examine the anomalous dimension matrix appropriate for the phase space restricted B — ► 
X u I v and B — > X s 7 decay spectra to subleading nonperturbative order. The time ordered prod- 
ucts of the HQET Lagrangian with the leading order shape function operator are calculated, as are 
the anomalous dimensions of subleading operators. We establish the renormalizability and closure 
of a subset of the non-local operator basis, a requirement for the establishment of factorization 
theorems at this order. Operator mixing is found between the operators which occur to subleading 
order, requiring the subleading operator basis be extended. We comment on the requirement for 
new shape functions to be introduced to characterize the matrix elements of these new operators, 
and the phenomenological consequences for extractions of \V u b\. 
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I. INTRODUCTION 



Extracting the CKM parameter \V u b\ is an important step in testing of the CKM de- 
scription of CP violation in the B meson system. Currently, the theoretically cleanest 
determinations of \V ub \ come from inclusive semileptonic decays which are not sensitive to 

V ub \, utilizing 



the details of hadronization; although recently an approach of extracting 
B — > 7T7T, has been advanced with a competitive error to inclusive methods. 

In inclusive extractions of \V u b\, experimental cuts to exclude the charm background of 
B X c £ v are imposed. This restricts the decay products to hadronic final states that have 



large energy E x ~ m B and low invariant ,„a S8 M x ~ ^ B A Q0D . With these phase spaee 
restrictions the local OPE expansion |2| appropriate for sufficiently inclusive decays used to 



extract \V r 



cb\ 



3, 



typically breaks down. 



As the local OPE, and the clean separation of scales that the local OPE represented in the 
analysis of B — > X c i v , is no longer valid; a more involved theoretical approach is required 
to separate the scales relevant to these decays. Decay rates are expressed as convolutions of 



hard (if), jet (J) and soft physics (S) associated with the scales m b ^> ^/Aqcd^6 ^ Aqcd, 
in the following way, 



dF = H^,a s (^ jduji^^M^ s {<*>)• (i) 

Although this factorization theorem has been proven diagrammatically j^] at leading order 
in 1/rrib, it is not known to hold to all orders in the nonperturbative expansion. It has only 
recently been extended beyond leading nonperturbative order |(|. 

The systematic treatment of the nonperturbative corrections involve a two step matching 
procedure. One matches QCD onto the effective field theory of the intermediate scale, 
describing quarks and gluons with large energy and small offshellness, known as SCET 
I23I 24, 25 1 , and uses the renormalization group evolution to run down to the soft scale. 
One then matches SCET onto the lightcone wavefunction of the B meson, expressed in 
terms of HQET fields. One can also match directly from QCD onto HQET, a much simpler 
procedure at the cost of not summing the logarithms of the ratio of scales log ^771^^^/771^ 
via SCET. In either case, the soft sector of the theory is expanded in terms of non-local 
operators. The leading order term in the AqcD/^b expansion of the lightcone distribution 



function of the B meson 
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is know as the shape function. 



At subleading order in the nonperturbative ex pan sion, four additional non-local operators 
have been determined to be present 0, 10, [ill I12L Q], the matrix elements of which are 
referred to as subleading shape functions. 

It is of some intrinsic interest to examine the renormalization of these non-local operators, 
as they are non-local and their renormalizability is not know a priori. It is also important 
to know if this set of operators is complete in the error assigned to extractions of \V u b\, 
and in considering the above factorization theorem beyond leading order. Examining the 
perturbative behavior of these subleading shape functions is also a necessary step to take 
to perform one loop matching calculations onto the soft sector. For these reasons we have 
determined the anomalous dimension to subleading nonperturbative order. 

To this end, we have determined the contributions of the time ordered products of the 
subleading £hqet with the leading order shape function, and examined the anomalous 
dimensions of the subleading nonperturbative operators. We establish the renormalizability 
and closure of a subset of the subleading non-local operators. We find that the known 
operator basis mixes with new operators, requiring that the subleading operator basis be 
extended. We also comment on the phenomenological consequences of these results. 

II. ANOMALOUS DIMENSION TO SUBLEADING ORDER 
A. Notation 

We introduce two light-like vectors n M and n^ 1 related to the velocity of the heavy quark 
by v — 5(71 + fi), and satisfying 

n 2 = n 2 = 0, v ■ n — v ■ n — 1, n ■ n — 2. (2) 

In the frame in which the b quark is at rest, these vectors are given by n M = (1,0,0, 1), 
W 1 = (1,0,0,-1) and = (1,0,0,0). The projection of an arbitrary four-vector a a onto 
the directions which are perpendicular to the lightcone is given by a" = g^ap, where 

= g" v - - [n»n v + n u n») . (3) 

We also define a perpendicular Levi-Civita tensor 

ef = e al3cTp v a n p . (4) 



We also use the projector P + = 1/2(1 + as well as the Dirac structure s v = P + 'j Vr y 5 P + , so 
that v ■ s — 0. 



1. Distribution Notation 



Rather than the usual definitions of the star distribution as given in Neubert and deFazio 



'log(x) 



x 



lim 
lim 

/3-*0 



9{x-(3) 
x 

9{x-(3) 
x 



+ S(x - (3) \og(x) 



log(x) H — S(x — (3) log 2 (x) 
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we utilize the alternate notation, equivalent to the /^-distribution's in 

1 
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b n (x) = lim 



(x - (3) log 



n+l, 



X] 



n + l 

This notation has a fairly easy correspondence to the usual star distribution notation 

1\ , s /log(x)\ 



<P' (x) 



X 



<P[(x) 



X 



A useful identity given by analytic continuation is 

6{x) 1 



x 



l+e 



: 5(x) + 0' o (x)-e0' 1 (x) + O(e 2 



(5) 



(6) 



(7) 



(8) 



This relationship is valid when integrated against arbitrary functions f{x), where / is not 
singular at the origin. In general we can write the recursion relation 



9{x) 



1 d 



9{x) 



x 



n — l+e 



for n > 2. 



(9) 



x n+e n _\J r6 fix 

Several other useful properties of this function are (for some positive constant a): 

x <p' (x) = 6(x) 
x<f>o(x) = 5(x) - <f)' (x) 

a^' Q (ax) = <f>' Q (x) + 6(x)log(a), (10) 



as well as 



'9{x) 



f(x) + f{0) log(a). 



(11) 



/o \ x 

For an a with dimension, such as a = Aqcd, the above equation is suitably modified to 
obtain logs the dimensionless ratio a//i. 



B. Operators to subleading Order 



At leading order a single non-local operator characterizes the nonperturbative physics, 

Q (u,T) = h v 5{u + m-D)Th v , (12) 



where the covariant derivative is = <9 M + ig s -^-n- 

The order Aqcd/^6 corrections to the B — > X u £u and B — » X s y decay spectra require 
the introduction of four additional non-local operators js, 10 



-» Xj dec 



m b Qi{ui,Y) = h v {iDj_, 8{uj + in-D)} Th v , 

m^Q^iyj.Y) = h v [iDj_,5(u + in-D)} T h v , (13) 
iTii, Qz{uj, r) = J dix)\dix)2 5{uj\,uj 2 ] u)h v 5{u 2 + in-D)gj"{iD±, iD^Siui + in-D)Y h v , 
nib Q<i(u, T) = —J duj\duj2 5{uJxi lu 2 ; u)h v 5(tu 2 + in- D)iej^[iDj_, iD±]5(uJi + in-D)Y h v , 
where 

5(u -wi) - 5(u -u 2 ) , s 

d(a;i, w 2 ; = • (14) 

LO\ — UJ 2 

We define these operators rescaled by m& for later convenience in the anomalous dimension. 
This rescaling should be noted when comparing to other work dealing with subleading shape 
functions. We also use the convention of labeling operators as Qi operators when the Dirac 
structure is general, and referring to the operators as 0« and when a particular Dirac 
structure is required. 

We find that the operator basis must be extended beyond tree level to include, at least 
the following operator 

m b Qi(uj,T) = —2l ' duJiduj 2 8(iUi,LU 2 ;u;)K 2 (u 1 ,u 2 ;T), (15) 

where we have defined the following kernel and coefficient functions 

K 2 (uji, u 2 ; r) = h v 8(uj\ + in- D)iD fJ L 5{oj 2 + in-D)Th v , (16) 

e^^u,) ^-^-^-^ . (17) 

tf - "* nn 

The operator Q\ was originally defined in p| by Mannel and Tackmann [2$, |2J| based 
on symmetry arguments and examining the endpoint of B — ► X c I v and taking the massless 
limit. We find that the operator is unambiguously required beyond tree level due to the 
mixing experienced with the original set of operators. 



C. Operator Feynman rules 



We use Feynman Gauge to calculate the anomalous dimension to subleading order as the 
usual choice of lightcone gaug e introduces non physical poles in the calculation, for a review 
of the relevant issues see |30j. 

Below, we present the required one and two gluon Feynman rules. The two gluon Feynman 
rules are also required and used but are too lengthy to include here, they can be obtained 
from the authors upon request. The non- vanishing zero gluon Feynman rules in Feynman 
gauge are: 

(h v (k)\Q (u,T)\h v (k)) = 5(cu + n-k)T, 
(h v (k)\Ql(u,T))\h v (k)) = 2^5(oj + n-k)T, 

{h v (k)\Ql(u,T))\h v (k)) = 2^5(oj + n-k)T, 

m b 

{h v (k)\Q 3 (uj,T))\h v (k)} = -2^6'(uj + n-k)T. (18) 

m b 

The one gluon Feynman rule for the leading order operator is 

(h v (k)A:(£)\Q (u;,T)\h v (k)) = -gT a n v V. (19) 

The one gluon Feynman rules for single covariant derivative operators are: 

(h v (k)A:(£)\Ql(u;,r)\h v (k)) = -gT a g^6 + (n-£)— - g T a n» (2k + £)» x ( S -^^) — , 

m b \ n-i J rrib 

<M*)*MI*<«.r>M*» = -w(^)^ + 2 9 r^(^)^ 

\ n-i J rrib \ n-£ J rrib 

(h v (k)A:(£)\Q%(cu,r)\h v (k)) = gT a g>?8^n-£)—-gT a n v e i _ — . (20) 

rrib \ n-i J rrib 

Finally, the one gluon Feynman rules for two covariant derivative operators are as follows: 
(h v (k)A:(£)\Q 3 (u,r)\h v (k)) = 2gTj(2k + £)l( 6 -^^)-n»k 2 f^^ + ^^W T 



(h v (k)A» a (£)\Q A (u,T)\h v (k)) = 2gT a ief£ ±fj — . (21) 
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where t is the gluon momentum flowing out of the vertex, and the gluon carries Lorentz 
index v and colour index a. We have also made the convenient definitions 



6±(x) = 9(u + n-k + x)±9(uj + n-k) (22) 
6±(x) = 5(uj + n-k + x) ± 5{uj + n-k). (23) 

D. The Anomalous Dimension Matrix 

The renormalization of the operators Qi(u, T) is performed in the usual fashion, 

<2i(^,r) barc = J du'Z ij (uj',uj,fl)Q i (u',fl, r) ren , (24) 

where Zij(u',u,fx) is a matrix of renormalization constants. The values of the elements 
of Zij can be found by taking arbitrary partonic matrix elements of both sides, which at 
leading order gives Z^\lu', lu, fx) = 5ij5(uj — u'). 
To subleading order in a s we have 

(g i (^,r)) b ° a ) rc + a s (Q i (^,r)) b 1 a ) rc = [du/\Z§\u\u,fi) + a a zg\u\u,fi] (25) 



<g i (o/,/i > r)>£ + a.<g i ( £ «;',/i,r)>g 



) 

ren 



from which one obtains 



du'Z%\u\u^)(QA"'r^)) m = (Q i (uj,T))^ rc -(Q l (uj,T))^ n 

= ((Qi(^,r)) b 1 a ) rc )div (26) 

where by ((Qi(u, r)) barc )div, we refer to the UV divergent part of (Qi(u},T))^ Te . Because 
there are operators such as Q2 and Q4 which do not have a zero gluon form, we must consider 
matrix elements of Eq. (j24j) with at least one external gluon. These will be sufficient to 
identify the mixing of the various operators into Q2 and Q4. It should be noted that matrix 
elements with zero and one external gluon states are not sufficient in principle to determine 
the anomalous dimension matrix to subleading order. The operator 

O^fa, ou 2 , T) = h v [iD*i, 5{uo 2 + in-D)] [iD u ± , Sfa + in-D)} T h v , (27) 

does not have a zero gluon or one gluon Feynman rule. Its nonvanishing Feynman rules start 
at two gluon external states. In this paper, we will not be calculating the two external gluon 
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diagrams necessary to find mixing into this operator, if any exists. We extract the anomalous 
dimension matrix of the subleading operators by examining matrix elements containing one 
perpendicularly polarized external gluon: 

The non-perpendicular components of the gluon field were also examined but found to induce 
no further mixing. 

The mixing of Qq into the other operators is determined by calculating matrix elements 
of this operator with insertions of the subleading HQET Lagrangian. Zero gluon matrix 
elements are sufficient to find the mixing into Qo, while one gluon matrix elements are 
required for mixing into the remaining operators. Due to the spin symmetry violating 
effects of the subleading HQET Lagrangian, the anomalous dimension of the Pj operators 
can differ from that of the Oj operators. 

The wavefunction renormalization of the bare operators expressed in terms of renormal- 
ized fields are Qi(u, r) barc = Z h Z^ 2 Qi(uj, T) where n is the number of gluons in the operator, 
and Qi(u, T) is written in terms of renormalized fields. For diagrams with an external state 



gluon we use the background field method to treat the external gluon as a classica i 
so we aquire no Z 3 factor due to the wavefunction renormalization of the gluon. 



field and 



3l| 



E. Diagram Calculations 

One Gluon Matrix Elements 



The one gluon matrix elements are determined by calculating the diagrams shown in 
Figure for each operator. The external gluon in each of these diagrams is a background 
field gluon and the external states are chosen to have perpendicular polarization. We utilize 
dimensional regularization and the MS scheme to regulate our divergences. To isolate and 
remove the IR divergences in the calculation we keep all the particles off shell by retaining 
factors of v-k, v-l and £ 2 , where I is the external gluon momentum. 

To clearly illustrate the need to extend the operator basis we present the results for Qi 
diagram by diagram. In general, for perpendicular polarized external gluons, only diagrams 
(ac), (ad), (be), (bd) and (dc) contribute to the amplitude. For diagram dc, the loop 
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(aa) (ac) (ad) (dc) 



3 



(bb) 



i 

(be) 



■ a 



(bd) 



1 



(dd) 




(ca) (cb) (cc) (cd) 



FIG. 1: The one gluon diagrams which must be calculated for each operator. 

integrals to perform are as follows, with c\ = (a s g s T a g^) /(47r), 
,(i) _ c i n e\2 .A-d f A S+(n-£) 



mi, J {27T) a n-q [n-t + n-q) (q z + ie){{q + iy + le) 



+ —C A (n-£) 2 ^- d f 
rrih J 



mb J (27i) d n-q (n-£ + n-q) (q 2 + ie)((q + £) 2 + ie) 

The integrals are perfomed via the standard techniques of dimensional regularization with 
d = 4 — 2e, the MS renormalization scale ft 2 = 4iru 2 e~' yE , and the utilization of Eq.(jHJ and 
we have suppressed the lorentz and colour indicies. The UV poles obtained for this diagram 
for Qi after the integrals are performed and consideration of the symmetry factor are 



(*Aic) 



(!) c\ Ca I n-i , u + n-k + n-£ n-i , u + n-k, 

div = m h e \p,{u + n-ky ^ p, ' ~ ~u(uj + n-k + n-ty^ Ji ' 



. c i c a ( x f «\ i ,n-i y 

+ 5 + {n-£) log(— , 

m b e \ u 



The results for diagrams ac and be when inserting Qi are 

... .m ci Ca f S(oj + n-k + n-i) 

(zAc)Si = o - + 5(uj + n-k + n- 

mbe \ 2e 
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ci C A (, ,n-i s 

log(— — ) o(lj + n-k + n-i 



m b e V /i 



CiCa / (cj + n-A;) uj + n-k uj + n-k + n-£ \ 

-v{ — = — = ) , (duj 



m b n-£e y (a; + n-A; + /i /i 

/• * an C!Ca f 5(lu + n-k) ,n-£. ,.\ 

?7ifc e y z e /i j 

c x C A ( (u + n-k + n-i) u + n-k + n-£ u + n-k\ 

H 7" 7 — ; n — y l = > + y l = ) ■ K 61 ) 

m b n-ie\ [00 + n-k) fi (x J 

Finally, the results for diagrams ad and bd for Q\ insertions are 

,. - v m ciCf (5+(n-i) 2 . uj + n- k + 2 ...u; + n-fc x \ 

Mad)div = — — TM ~ ) _ ~ Po( 



+ 



m^eye /i /x /i /j 

c x C A ( 5(u + n-k + n-£) (f)' ((uj + n-k + n-i)/^ 
mje I 2 e /7 , 



,. . v m ciCp ( 6 + (n-£) 2 a> + n- fc + n-f. 2 a> + n-fc. A 

«Ad div = -0o = - -M — = — 

m b e y e /j /i /i A* / 

ciCU / 5(uj + n-k) <p' ((uj + n-k)/(l) \ 
m b e y 2e /I y 

The amplitudes combine to give the following UV poles 

„ = W _ A^n- m \ + / + 29. { n-m\ 

mv m b e \ e jl J m b e \ n-£ J 

Once the wavefunction renormalization terms are multiplicatively combined with the result, 
we express the amplitude in terms of renormalization matrix elements dj and the one gluon 
Feynman rules for the operators Qi and Qi as follows 

m«i>£ = j-J Md x {u,tj,m<%w) m (34) 

a 



+ ^ l duj'd,{uj,uj') (<0f(o/)> (0) - (Q5V)> (0) ) 



where 



, / / ~\ 2(7p . 2C^ . .. AC p . ( uj — u/\ 

di{uj,uj ,fjL) = — d[U) — v)-\ d(uj-uj) + — — O — - — , 

e z e fie \ jl J 

d 4 (uj,uj') = —^-8{uj — uj'). (35) 

The form of the mixing between Qi and Qi deserves some comment. At zero gluon the 
matrix elements of these operators are identical causing this mixing to be undetermined for 
zero gluon external state diagrams, even though the zero gluon matrix elements of both 
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operators are nonzero, contrary to naive expectations. The contribution of the operator Q± 
to the renormalization matrix was also determined. We find that this operator mixes with 
itself contributing a d\ form to the matrix Z$l- The antisymmetric operators Q2 and Q4 mix 
only with themselves and contribute diagonal factors of d\ to the matrix of renormalization 
constants. 

The inauspicious form of mixing between Q± and Qi is also present in the operator Q 3 , 
however the corresponding analysis of Q3 is more complicated and is still under investigation. 
Due to this complication in determining the full anomalous dimension matrix and the need 
for a two gluon calculation to determine the possible mixing with Q^^ui, U2, T) we present 
the results of our initial study of the anomalous dimension to subleading order in this paper 
and comment on the phenomenological consequences of the presented results. We collect 
our results in section IIHI 

The T products of Chqet with Oq and Pq 
To find the mixing of the operators 

O = h v 8{uj + in- D) h v , 

P M = h v 8(u + in-D)j^ 5 h v (36) 

into the subleading operators, we must evaluate the time ordered products of the operators 
with the the subleading terms of the HQET Lagrangian (Ci) 

T {uo) = 1 J d i xT[C 1 (x),O (0,uj)}, 

T£(uj) = i y'd 4 xT[£ 1 (x),P M (0,cu)]. (37) 

We now explicitly refer to the Dirac structure of the operators. This is necessary due to 
the Dirac structure of the operators in the subleading Lagrangian. We treat the subleading 
Lagrangian as a single operator insertion for the purposes of our calculation. The different 
renormalization of the kinetic and chromomagnetic terms is accommodated by breaking the 
T products up in to 7(o 0) o fc )> ^(Oo,o m )> T(p 0t o k ), T(F ,o m ) after the diagram calculations, where 
Ok, O m refer to the kinetic and chromomagnetic operators of the subleading Lagrangian. 
We start with the zero gluon diagrams. They are illustrated in Figure El The crosses in 
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On 



On 



On 





(a) 



(b) 




FIG. 2: The diagrams which must be calculated for Qq. The crosses represent insertions of the 
subleading HQET Lagrangian. 



the diagrams represent the locations where one inserts the subleading HQET Lagrangian, 
given by 

_ (n n , \ 2 _ (nrr „ H a P\ 

(38) 



>"i = v ^ h>v ~ cL{fl)n v h v . 



2mf, Arrib 
The zero, one and two gluon Feynman rules for this Lagrangian, where we suppress the 

renormalization scale dependence of the O m operator, are 



z£i[0-gluon] 
i£L\ [1-gluon] 



2rrih 



-igT a 



(2k ± + 



2m>h 



-P+ + igT a - 



2mh 



OL\CC2 

iA[2-gluon] = ig 2 {T a \T a *} 9 -^P + + tg 2 [T a \T a2 ]- 



:a 1 a 2 pri 



(39) 



2rrib 2rrib 
The UV divergent part of the sum of the subleading £hqet zero gluon results for both 
Oo and Pq (up to Dirac structure) is the same. Our result for the general Dirac structure 
operator Q , with indices suppressed, is 

(i) a s ^ f j r Au'5(uj - u') 



(hy{k)\T {QoPk) {U,T)\hy{k))^ y 



ArribTT 



C F / du 



■(hy(k)\Q (cU>,T)\hy(k))W 



+ ^L Cf f du' 36{U ^ v^hyik^Q^uo'J^Kik))^. (40) 
47r J e 

The latter term is the zero gluon matrix element of a modified Q± operator. The operator is 
modified to not have a perpendicular covariant derivative, but simply a in its definition. 
This term vanishes at this order due to the equations of motion, these results are combined 
with our operator anomalous dimensions in section 11111 (It should be noted that the zero 
gluon calculation does not determine this mixing is with the Q\ operator, as its zero gluon 
rule is identical to Q\. Consistency between zero and one gluon external state calculations 
determines this mixing to be with the modified Q\ operator.) 
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There are many more one gluon diagrams than zero gluon diagrams, as illustrated in 
Figure H3 The diagrams explicitly given constitute half of the total number of diagrams that 



> X X g X ! > > g X ■ > > X X =^ 



I 

(aa) (ac) (ad) 



ejSOTfc, r ^%, 



=*= > XXXX B ; > > XX B ^ =^= 



i i i 

c c °) 



(ae) (ca) (cb) 



if r 



(cc) (cd) (ce) 



FIG. 3: The one gluon diagrams which must be calculated for Oq and P$. The crosses represent 
insertions of the subleading HQET Lagrangian. The mirror diagram corresponding to each of the 
above diagrams is not shown. 

must be calculated for each of Oq and Pq. The other diagrams can be looked upon as either 
the mirror diagrams of those given, or the transposed diagrams which have the C\ operator 
and lightcone operator interchanged. 

The Dirac structure of the subleading Lagrangian force us to consider the Dirac structure 
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of these diagrams. Let us consider the one gluon diagrams of Figure El where the lightcone 
operator is Qo- We will denote by (..|Ar|--) and the amplitude of these diagrams 

and the amplitude of the mirror diagrams respectively. Because of the simple relationship 
between Oq and P£, the corresponding amplitues for Pq are (..|s°«4r|..) and (..\ALS a \..). 

Each of (A R ) and (Al) can be decomposed into the two Dirac structures P + and s v , for 
example with a heavy quark target, with one gluon in the final external state: 

(h v (k)A»\A R \h v (k)} = (h v (k)A»\A R+ P + \h v (k)) + (h v (k)AZ\Al a s r ,\h v (k)), 

(h v (k) A v a \A L \h v {k)) = (h v (k) A:\A L+ P + \h v (k)) + (h v (k) A:\Al s s v \h v (k)). (41) 

Thus for operator O we can write the total amplitude proportional to each of the Dirac 
structures after the calculations of the 40 diagrams required. The results for insertions of 
the Ok are: 

(K(k) A:\T i0 „ 0l )\K(k))2 = {K(k)A:\(A R+ + A L+ )P + \K(k))2, 



(42) 



nrribe \ n-£ 

3 C F a s9s ^ / M) _ n „ (2 k v + £ v) / 5 ± ^)\ \ p+ 



47rmb€ \ \ n-i 

This result is easily matched, it is identical to the mixing form found in the zero gluon result 
with r = P + , so that the mixing occurs with the operators 0$ and 0\. 

(hv(k) A:\T {O0>Ok) \h v (k))2 = ^C F [ du' 4 ^ ~ U ' ] (h v (k) A:\O (uj')\h v (k)) m 

+?C F / dJ 3 6{U - U>) v, ( K ( k) A: I Of (a/) \h v (k)) W . (43) 
47r J e 

The result of the T product with O m is 

(K(k) ^|r (O0lOm ,|A,(*)>£ = (*,(*) A" a \\Al, + Al,]s,\K(k))2, 

- C ^ aT '< i ^_(n. t )-i^ S -^f\. (44) 



8 nib en y n-i 
This result matches onto the one gluon rule for P 2 , as expected by the symmetry of the 
single derivative operators: 

(K(k) A:\T {O0 ,o m) \h v (k))2 = -%-C A J d J^'^ ie^ (h v (k) A:\P 2 ^>)\h v (k))W . 

The total T° Pq q ^ and T? Po 0m j amplitudes can be written as 

(h v (k) A:\Tfa h0 Jh v (k)) = (h v (k) Al\s°A R+ \h v {k)) + (h v (k) A» a \s°A L+ \h v (k)), (45) 
(h v (k)A:\Tfo >0 Jh v {k)) = (h v (k)A:\s° S i(A Rs ) v \h v (k)) + (h v (k) A v a \s n s a {A Ls )r)\K{k)). 
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Using the decomposition 



s s 



(46) 



we can decompose in terms of the pieces proportional to P + and s v for these amplitudes. 
For Tfp a \ and T° Pq q \ we find the following mixing 



(h v (k) A:\Tfr 0>Ok) \h v (k))V = (h v (k)A:\s°(A R+ + A L+ )\h v (k)) 

,4u'5(u - uj') 



= —C F / duj 

4 71 
4tt 



em b 
,3 5{uj — uj') 



(h v (k)A:\p (ur\K(k)) 



v^(h v (k) Al\P^{u'Y\h v {k))^\ 



(hv(k) Al\T^ Jh v {k)) 



(i) 

div 



_ g^)(h v (k) A:\((A m ), + (A Ls ),)P + \h v (k)) 

+(h v (k) AZ\[i^Vr({A R8 ) n -(A La ) v )s < i > \h v {k)), 

^C A [ du' 6{uJ ~ ^ zer (K(k) A20 2 ,{u')\h v {k)) 
4 7r J Z e 

~ C A J duo' 6 ^~^ [{v° - rf)g^ + n V - w ff u")] x 



(MfcMZIi¥V)M*)> (0) - M*M^V)IM*)> 



III. RESULTS 



A. Leading Nonperturbative Order 



The order a s perturbative and leading order non-perturbative anomalous dimension ma- 
trix has been calculated by a variety of authors (itI 19^ . Our results agree with theirs, and 



in the basis 



the pertrubative expansion is given by 



{0 ,Po} 



(47) 



Z {1) (uj,uj',jl) 



5(uj - uj') 
5(uj - uj') 



Air 



di(uj,uj',n) 

di(oj,u',fl) 



(48) 
(49) 
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The distribution di(u, a/, p) is the combination of the operator and wavefunction renormal- 
ization counter terms, given by 



2C f 



2C 



5(u - J) + —5(u -cu') + - 



AC F , ( u -u r 



e l e fie 



(50) 



Recall, our initial expression related the bare and renormalized operators, 



Qo(u, r) barc = J dJZ{J , uj, p)Q Q {u\ f~L, T) 



(51) 



We differentiate this equation with respect to jl to obtain our renormalization group equation 



duj'j(uj',uj,fl)Q (uj',fl,r) ieri . (52) 

The anomalous dimension matrix is determined using the useful result for MS j^] 

dZi(a s ) 



■2a, 



da s 



(53) 



where Z\ is the coefficient of the 1/e poles. We find 

7(w,w',/2) = 

with, 



7i(u;,u/,/j) 

71 (w, a/,/*) 



7i(a;,a; ,/x) = CV 



7T 



c/ A 2 .. (UJ — Uj' 

5(u - u') + — 



(54) 



(55) 



B. Subleading Nonperturbative Order 

We have determined the matrix of renormalization constants at subleading non- 
perturbative order Z$l, excluding operators of class Q 3 . If we order our Qi operators 
as 



Oi = {O ,T(p ,o k ),T(p ,o m ),O%,O%,O!j;,O4}, 

T)<y J per rpa rpa p"f pffC po> r>a\ 

' i — VO > 1 (Po.Ofe)' 1 (fb.Om)' ^1 ' ^1 ' ^2 ) ^4 Si 



(56) 



the leading order term in the perturbative expansion in the basis {Oi,Vi} is given in block 
form as 



r Oj)O .(u;,a/) T° 0i>Pj (u,u') 

r° O >,o/) r° p>,u/) 



(57) 



1(3 



Where the entries in the matrices in the above expression with i,j = 0, 1...6 are given 



r p 4 A (<*>, a;') 



0. 
0. 



While the 0(a s ) term in the expansion is 



with the diagonal block matrices 



4tt 





d 2 (u;,u}',fl) d$(u,L/) 



d 4 (u},u',fi) d 5 (u,u') 

d x (u,u\p) 

di(w,o;',/i) 

di((j,u',p) 





di{uj , uj' , p) 










d£(w,w') 

d$ a71 {u,uj') 

d±{uj , lo' , p) d 5 (u,uj') 



di(u,u',fi) 












di(u;,a/,/i) 
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The off diagonal block matrices follows 



v o it pA u ,u') 







-d% a (u,uf) 











r ko,(^W) 







d% a {u,u') 

















The di(uj , uj' , jl) distributions are given by 

4CV u' 



d 2 (u,u ,p) 
d%(u,cv') 



e m b (p) 
3C F 



— 8{u — u>'), 



5{u - u/)if, 
C A 



d±(uj,uj',p) = di(u>, cu', p) 8{uj—uj'), 

d 5 (cu,uj f ) = —^-S(u — uj'), 



-^{ieT)5{u>-u>) 



Ca 
2e 



{{v a - n a )g^ + n V - v a v^))5(lu - u'). 



(60) 



We directly determine the anomalous dimension matrix to subleading order using Eq. ()53|) 
to be the following 

/ 



7o i! o,(^,^ / ,At) Tcvp^VjAi) 



(61) 
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The diagonal entries of the anomalous dimension matrix are 



j 2 (v,uj',fi) jg(uj,uj',ji) 



74(0;, cu', ji) 75(0;, cu', ji) 

7!(u;,u/,/I) 

71 (cu,cu',p) 





70iA (^^aO 










jx(oj,u',p) 





j 2 (oj,oj',P) 7^(0;, cu', p) 
















- 7 r>V,/i) 7? ff >X,/i) 
74(0;, a/, A) 75(0;, cu' , p) 









The off diagonal entries follows 



1fo i ,pAw,u',ii) 























ji(u,cu',p) 



c 
c 









71 (w, 











7i(w,a/,/2) 













lp,0, (v,u)', p) 







^ a {cu,cu') 
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With the following entries in the anomalous dimension matrix, 




72 



-3a s (p)—^5(uj - u') 

2 71 

, I ~\ , OL s [p) C A . , 

, n) + — 5{u -u), 

l 71 




74 (w, U)' 



P) 



7 5 (o;,u/ 




— - — 5{u-w), 

A TT 




77 K w 



A) 



(62) 



IV. CONCLUSIONS 

We have examined the anomalous dimension matrix appropriate for the phase space 
restricted B — > X u I v and B — » X s ^ decay spectra to subleading nonperturbative order. 
The effects of the time ordered products of the HQET Lagrangian with the leading order 
shape function operator were determined and the renormalizability and closure of a subset 
of the non-local operator basis used to describe these spectra, to subleading order, was 
established. 

Operator mixing was found between the operators which occur to subleading order, re- 
quiring that the subleading operator basis be extended to include the operator Qi- This 
requires the introduction of new shape functions to characterize the decay spectra of 
B — > X u i v and B — ► X s 7 beyond tree level. The mixing determined between the operators 
Qi and Q\ is of the pernicious form that required a one gluon external state calculation to 
determine, despite the non vanishing zero gluon Feynman rules of the operators. We have 
also demonstrated that the possible mixing with the operator Q^ioui, 002, T) in a similar 
manner; with vanishing Feynman rules for zero and one gluon, requires a two gluon exter- 
nal state calculation to completely determine the anomalous dimension at subleading non 
perturbative order. 

Mixing was also determined between the T product T(o .o k ) an d the leading order shape 
function, and the T product T(o ,o m ) was shown to lead to mixing between the Pj and Oj 
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operators at this order. 

The anomalous dimension and running of the Q±, Q% and Q4 operators was shown to be 
identical to the leading order shape function Q . 

This work can be built upon in a number of ways. The anomalous dimension of the oper- 
ator Q 3 is under investigation by the authors to establish the closure at one loop of the set 
of subleading operators discussed in this paper. The anomalous dimension of the sublead- 
ing four quark operators should be investigated to determine the full anomalous dimension 
matrix at subleading order. Once the full anomalous dimension is determined, Sudakov 
logarithms in the perturbative corrections to the subleading operators can be resummed, 
so that renormalization group improved calculations can be undertaken for the B — > X u i v 
and B — > X s 7 decay spectra to subleading nonperturbative order. 

More important than the small effect that these corrections have directly on the extrac- 
tion of \V u b\, is the fact that this work establishes the renormalizability of a subset of the 
soft sector nonperturbative expansion beyond leading order. This is a necessary step in ex- 
tending QCD factorization theorems beyond leading nonperturbative order, validating the 
factorization based approach used for the phase space restricted B — > X u ivi and B — > X s 7 
decay spectra beyond leading nonperturbative order. 
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